Better synchronizability predicted by a new coupling method 
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In this paper, inspired by the idea that different nodes should play different roles in network 
synchronization, we bring forward a coupling method where the coupling strength of each node 
depends on its neighbors' degrees. Compared with the uniform coupled method and the recently 
proposed Motter-Zhou-Kurths method, the synchronizability of scale-free networks can be remark- 
ably enhanced by using the present coupled method, and the highest network synchronizability is 
achieved at /3 = 1 which is similar to a method introduced in [AIP Conf. Proc. 776, 201 (2005)]. 
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I. INTRODUCTION 

Many collective dynamics in social, biological and com- 
munication systems can be properly described by com- 
plex networks. These networks exhibit complex topolog- 
ical properties such as the small-world effects and the 
scale- free properties Q, |^ H, Q . Many kind of network 
models have been made to embody these properties. The 
so-called small-world networks are the intermediates of 
regular lattices and random networks in structure but 
bear both characters of the two kind of networks, that 
is, they have small average distance as random networks 
and large clustering coefficient as regular ones 0. The 
scale-free networks are a kind of small-world networks 
with degree distribution obeying a power-law form. A 
scale-free network can be created by successively adding 
new nodes to the network and connecting them with the 
already existing ones by the preferential attachment rule 

H 

The interesting topological properties of complex net- 
works make the dynamics taking place on them much 
different from those on regular or random ones. For 
example, coupled dynamical oscillators on small-world 
networks are much easier to synchronize than on regu- 
lar lattices, and increasing the proportion of shortcuts of 
networks will make the oscillators more synchronizable 
0, 0, 0j [lO]- It bas also been observed that the more 
heterogeneous of the network degree distribution is the 
harder for the oscillators on the network to synchronize 
Therefore, generally speaking, networks with short 
average distance and homogeneous degree distribution 
will make the oscillators on them more synchronizable 

[Ulillillll. 

Very recently, motivated by practical requirements and 
theoretical interest, numbers of researches have begun to 
study how to enhance the network synchroniza bility. es - 
pecially for scale-free networks |T1 111 111 [ll |20| . 
The method proposed by Zhao et al. [l^ can sharply 
reduce the maximal betweenness thus enhance the net- 



work synchronizability, but it will bring some economic 
and technologic problems since the network structure is 
slightly changed. The method proposed by Chavez et 
al. keeps the network topologyunchanged, while adding 
some weight into the system 0,123. However, to com- 
pute the weight, this method needs the global structural 
information, which is usually unavailable in huge com- 
munication systems. Therefore, in this paper, we keep 
the network topology unchanged, and concentrate on the 
coupling method using only the local information. The 
Motter-Zhou-Kurths (MZK) method 15] is a typical ex- 
ample, in which the coupling strength from a node i is 
inverse to its degree fc^. In MZK method, every neighbor 
of a node has the same influence (coupling strength) to 
this node. However, in real networks, different nodes may 
have different influences. For example, in society, some 
people have strong influence on others in some aspect 
but they are not influenced at the same level. Another 
impressing phenomenon is that in the World Trade Web, 
the small countries' economies fluctuate with the power- 
ful countries tightly, but the contrary does not occur 2i|. 
Thus here, based on the assumption that different nodes 
play different roles, we adjust the influencing strength 
of each node receiving from their neighbors according to 
the neighbors' degrees. That is, a node is not influenced 
by its neighbors equally. It is found that oscillators on 
scale-free networks coupled in this way can have much 
stronger propensity for synchronization than in the pre- 
vious ways, with the exception of a recently introduced 
method [la. , where the former is similar to 
beta = 1 in our model and is shown here to be the most 
efficient method based on the information provided by 
the degree of nearest neighbors only. 
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This paper is organized as follow: in section 2, the dy- 
namical equations of coupled oscillators and the master 
stability function will be briefly introduced. In section 3, 
we will give the simulation and analysis about synchro- 
nization of correlated scale- free networks. Finally, we will 
draw our conclusion in section 4. 
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FIG. 1; (color online) (a) R in the parameter plane {a, (5). (b) i? vs a for different parameter /3. The numerical simulations 
are implemented based on the BA network of size A'^ = 1024 and with average degree k = Q. The data are obtained over 10 
independent realizations. 



II. THE DYNAMICAL EQUATIONS AND 
MASTER STABILITY FUNCTION 

For a network of N linear coupled identical oscillators, 
the dynamical equation of each oscillator can be written 
as 

i^ = F(x')-a^G,,H(x^"), z = l,2,...,iV, (1) 

where — F(x*) governs the dynamics of individual 
oscillator, H(x^) the output function, cr the coupling 
strength, and Gy the elements of the N x N coupling 
matrix. To guarantee the synchronization manifold an 
invariant manifold, the matrix G should has zero row- 
sum. Traditionally, the oscillators are coupled symmet- 
rically with uniform coupling strength and the coupling 
matrix G has the same form as Laplacian matrix L, that 
is, Gij — Lij, where 

{ki for i = j 
-1 forjeA, (2) 
otherwise, 

where ki is the degree of node i and is the set of 
i's neighbors. Because of the symmetry and the positive 
semidefinite of L, all its eigenvalues are nonnegative reals 
and the smallest eigenvalue Ao is always zero, for the rows 
of L have zero sum. And if the network is connected, 
there is only one zero eigenvalue. Thus, the eigenvalues 
can be ranked as Aq < Ai < A2 < ... < Aat-i. According 
to the criteria of master stability function [2^ |2^ I25ll2^ . 
the network synchronizability can be measured by the 



eigenratio R = Xm-i/Xi'- The smaller it is the better the 
network synchronizability and vice versa. 

It is later found that networks with high heterogeneity 
of degree distribution coupled uniformly are hard to syn- 
chronize. As mentioned above, to eliminate this problem, 
Motter, Zhou and Kurths suggested the coupling matrix 
taking the form dj = U.j/k^ El El El. This simple 
change of coupling matrix enhances the network synchro- 
nizability sharply, and the optimal condition is = 1. 
And by exploiting the information contained in the load 
of each edge (i.e. set the off-diagonal elements of the zero 
row-sum coupling matrix G to be Gij = Ifjl ^j^ilfj, 
where lij is the load of edge connecting node i an d 7), fur- 
ther enhancement in synchronization is achieved |l9ll20j |. 

Many real-world networks are highly heterogeneous 
with a few nodes, named hubs, having very large de- 
grees. When using the uniform coupling method, these 
hubs synchronize first, and slowly the nodes with fewer 
degree synchronize to them 0| . If the influence of the 
hubs on the low-degree nodes becomes stronger, the lat- 
ter will synchronize to the former much easier, obviously, 
the network synchronizability will be enhanced. There- 
fore, we argue that not only reducing the communication 
load of hubs (as did in MZK method) , but also increase 
their influences may further enhance the network syn- 
chronization. 

Here we take into account the effects of different de- 
grees of nodes on synchronization, that is, a node in com- 
plex network is not coupled uniformly by its neighbors 
but the coupling strength is modulated by /c". Thus the 
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FIG. 2: The sketch maps of three simple equivalent networks, where the arrow form node i to node j indicates the latter 
receives coupling signal from the former. Their eigenratios are 2 (a), 6.8284 (b) and +oo (c), respectively. 

where D = di&g{k:["S:['^ ^k^'^S^^ , ...,kJj"S~'^} is a di- 
agonal matrix, and L' = {L[j) is a symmetric zero row- 
sum matrix, whose off-diagonal elements are L'^j = kfk'^. 
From the identity 

Qlci{DL' -\I)=det{D^L'D^^ (5) 

valid for any A, where "det" denotes the determinant and 
/ is the N X N identity matrix, we have that the spectrum 
of eigenvalues of matrix G is equal to the spectrum of a 
symmetric matrix defined as 

H = DH'Di. (6) 

As a result, the eigenvalues of matrix G are all nonnega- 
tive real and the smallest eigenvalue is always zero. 



III. SIMULATIONS 




FIG. 3: (color online) The eigenratio R vs parameter a at 
/3 = 1.0 for several BA network configurations of size — 
1024 and with average degree k — 6. Each color represents 
one configuration. 



coupling matrix G takes the form 

r SJS^ forz=j 
G,, = <^ -kpS^ for J e A, (3) 
[ otherwise, 

where Si = J^jeAi ^f- When a — P — 0, this coupHng 
scheme degenerates to the uniform coupling scheme , 
the case of a = corresponds to the MZK method |l5j |. 
and the case of /? = 1 is equivalent to the one introduced 
in the ref. in (see the Eq.(15) for details). 

Using a similar method to the one proposed by Mot- 
ter et al. we next prove that all the eigenvalues of 

matrix G are real. Note that, Eq. (3) can be written as 

G = DL', (4) 



In our coupling method, giving the parameter /3, for 
a > 0, nodes with large degrees have stronger influence, 
and for a < 0, nodes that bear few edges are more influ- 
ential. Parameter f3 is exploited to eliminate the discrep- 
ancies between the coupling signals that each node re- 
ceive: Given a, when (3—1, each node receives the equal 
quantum of signals, when /3 < 1, nodes that have larger 
sum of neighbors' degrees are influenced more strongly, 
and when /3 > 1, the contrary situation occurs. 

Figure 1(a) shows the numerical values of eigenratio 
R on the parameter space {a, 3) for the well-known 
Barabasi- Albert (BA) networks [g- To clearly exhibit 
the effects of a and /3 on R, we report the values of R as 
a function of a for different f3 in figure 1(b). No matter 
what value the parameter (3 takes, there exists a region 
of a, in which the eigenratio R is smaller than that of 
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FIG. 4: (color online) The eigenratio R vs the parameter a 
at l3 — 1.0 for the generalized BA networks with different 
assortative coefficients r. In all cases, the average degree is 
k = 6, and the network size is A'^ = 1024. The data are 
obtained over 10 realizations of network configurations. 
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FIG. 5: (color online) Ratio R vs the parameter a at /3 = 
1.0 for the generalized BA networks with different assortative 
coefficients r. In each plot, the average degree is fc = 6, and 
the network size is A = 128 (a), A = 256 (b), A = 512 (b), 
A — 2048 (d). The data are obtained over 10 realizations of 
network configurations. 



FIG. 6: The crossed value of a vs the network size A at 
l3 = 1.0 for the generalized BA networks. In all cases, the 
average degree is k — 6. 



the disconnected case, the eigenratio R will approach 
to infinite, while in the connected case, the eigenratio 
equals to 2 or some other larger constants. Figure 2 illu- 
minates three simple equivalent networks with a = +oo; 
the former two are connected, and the third one is dis- 
connected. Their eigenratio are 2, 6.8284 and +oo, re- 
spectively. Figure 3 shows the changes of eigenratio R 
with the parameter a with f3 — 1 of different network 
configurations. When a > 0.4, the eigenratios for differ- 
ent configurations go apart: Some approach 2 or other 
constants not much larger than 2, while some go to in- 
finity, due to whether the equivalent networks are con- 
nected or not. In addition, from the simulation, we find 
that with the increasing of network size, the proportion 
of networks being disconnected when a — oo (— oo) will 
increase sharply. 

Next, we investigate the effects of degree-degree cor- 
relation on the network synchronizability (29ll] . The cor- 
related networks are generated by an extended BA algo- 
rithm 30, 31]: Starting from mo fully connected nodes, 
then, at each time step, a new node is added to the net- 
work and m (< mg) previously existing nodes are chosen 
to be connected to it with probability 



P i oc 



the case a = 0. That is to say, when proper parame- 
ters are chosen, our coupling method can be even better 
than the MZK method. Similar to the results obtained 
from MZK method, f3 = 1.0 corresponds to the optimal 
case (i.e. the highest synchronizability). Hereinafter, we 
concentrate on the case of (3 — 1. 

Note that, in the limit a — +oo (— oo), each node is 
only influenced by the neighbor having the largest (small- 
est) degree. The similar situation as mentioned in Ref. 
[19j appears: The original network approaches to a new 
configuration that is connected by some effective directed 
edges p^ . and the new network, named the equivalent 
network, may be either connected or disconnected. In 



(7) 



where pi and ki denote the choosing probability and de- 
gree of node i, respectively. By varying the free param- 
eter ko (> — m), one can obtain the scale-free networks 
with different assortative coefficients r ^2t iMj • 

Fig. 4 shows the relationship between eigenratio R and 
the parameter a for different assortative coefficients given 
/3 — 1. Interestingly, there exists a unique cross point at 
aero ~ —0.25. When a < aero, the stronger assortative 
of network predicts better synchronizability, while when 
a > aero, contrary phenomenon appears. In Fig. 5, we 
report the simulation results for networks with different 
sizes, which heightens the reliability of the existence of 
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this crossed behavior. Fig. 6 exhibits the cross point aero 
as a function of the network size N: It increases when N 
is small, and will get steadily at about 0.25 for sufficiently 
large N. Although it is interesting, unfortunately, we are 
not able to provide a theoretical explanation about this 
phenomenon. 

IV. CONCLUSION AND DISCUSSION 

The hub nodes of a highly heterogeneous network al- 
ways play the major roles in determining the dynamical 
behaviors of the network. In synchronizing process, the 
hub nodes simultaneously have two effects. On the one 
hand, the throughput of these hub nodes are too heavy, 
thus they will hinder the coupling signals' transmission. 
On the other hand, they have great controlling capability 
for their large number of coupling neighbors. The MZK 
coupling method has taken into account the former point, 
and can predict much better synchronizability than the 
uniform coupling method. The present method further 
considers both aspects, performs even better than MZK 
method, and shows that synchronizability is maximum 
for a set of parameters that is equivalent to the method 
introduced in Ref. 0|. 



Some previous works [3 HE HE 113 suggested that 
there exists some essentially common features between 
network traffic and synchronization on a dynamical level 
since the performance of them both are mainly deter- 
mined by the maximal betweenness, and many methods 
used to enhance the network synchronizability can also 
improve the traffic conditions. Therefore, a natural ques- 
tion raises: will the network throughput increase if each 
node tends to receive information packets from the large- 
degree neighbors? 
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